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Abstract. We study perfect fluid cosmological models with a constant equa- 
tion of state parameter 7 in which there are two naturally defined time-like 
congruences, a geometrically defined geodesic congruence and a non-geodesic 
fluid congruence. We establish an appropriate set of boost formulae relat- 
ing the physical variables, and consequently the observed quantities, in the 
two frames. We study expanding spatially homogeneous tilted perfect fluid 
models, with an emphasis on future evolution with extreme tilt. We show 
that for ultra-radiative equations of state (i.e., 7 > 4/3), generically the tilt 
becomes extreme at late times and the fluid observers will reach infinite ex- 
pansion within a finite proper time and experience a singularity similar to 
that of the big rip. In addition, we show that for sub-radiative equations of 
state (i.e., 7 < 4/3), the tilt can become extreme at late times and give rise 
to an effective quintessential equation of state. To establish the connection 
with phantom cosmology and quintessence, we calculate the effective equation 
of state in the models under consideration and we determine the future as- 
ymptotic behaviour of the tilting models in the fluid frame variables using the 
boost formulae. We also discuss spatially inhomogeneous models and tilting 
spatially homogeneous models with a cosmological constant. 



1. Introduction 

In cosmology it is essential to specify a set of observers, or rather a congruence 
of world-lines, from which observations are made. Physical quantities, such as the 
Hubble parameter, depend on the choice of congruence and consequently the Hub- 
ble parameter is observer-dependent. In this paper, we will consider a universe with 
two different congruences of observers, with respect to which the interpretation of 
the Universe can be radically different. Spatially homogeneous perfect fluid uni- 
verses in which the fluid is not necessarily comoving with the spatially homogeneous 
hypersurfaces are prime examples of this phenomenon, and have the advantage that 
they are dynamical solutions to Einstein's field equations and can model some fea- 
tures of our real Universe. For example, in these models the observers moving with 
the geometric congruence can experience an ever-expanding non-inflating universe, 
while observers moving with the matter congruence can experience an inflating 
universe. 

For spatially homogeneous (SH) Bianchi cosmological models, the universe is 
foliated into space-like hypersurfaces (defined by the group orbits of the respective 
model) ^ El El • For these perfect fluid models there are two naturally defined 
time-like vector fields (i.e., congruences): the unit vector field, n a , normal to the 
group orbits and hence orthogonal to the surfaces of transitivity (the 'geometric' 
congruence), and the four- velocity, u a , of the fluid (the 'matter' congruence); see 
Figure 1. If u a is not aligned with n a , the model is called tilted (and non-tilted 
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Figure 1. The normal and the fluid congruences defined by the 
vector fluids n a and u a , respectively. The space-like hypersurfaces 
£ are orthogonal to n a and are defined for each Bianchi types by 
the group orbits of the respective Lie group. 

or orthogonal otherwise) |5J- The geometric congruence is necessarily geodesic, 
vorticity-free and acceleration-free. The matter congruence , on the other hand, is 
not necessarily geodesic and can have both vorticity and acceleration. 1 Usually, 
the kinematical quantities associated with the normal congruence n a of the spatial 
symmetry surfaces, rather than the fluid flow u a , are used as variables (a comoving 
fluid description is presented in Appendix A). Following [5J, tilt variables v a are 
introduced, so that in an orthonormal frame where n a = (1, 0, 0, 0), we have 



where v 2 = v a v a ,v < 1. 

We will assume a perfect-fluid matter source with p = (7 — as equation 
of state, where /1 is the energy density, p is the pressure, and 7 is a constant. 
Causality then requires 7 to be in the interval < 7 < 2. A positive cosmological 
constant may also be included in the models. Such SH tilted cosmologies with a 
7-law perfect fluid source have been studied by a number of authors. The tilted 
Bianchi type II model was studied in [S] , the type V model in [HI ED E3 \Ui 1 type 
VI in QjOUHlHl], types IV and Vll h in [151 ITS] , a subset of the type VI h models 
in ^3 , irrotational type VIIo models in El E] > general tilted type VIIo models in 
[T5| . and tilted type VIII models in [T5| . 

It is known that the tilt can become extreme asymptotically (i.e., v 2 — > 1). 
Although it is possible for tilt to be extreme both to the past and to the future, 
we shall focus on the late-time behaviour of the tilted Bianchi models here. The 
extreme tilt future attractors in SH models (without a cosmological constant) are 
summarized in Table 1 (using the notation of El ~ El)- I n particular, let T 
be the proper time as measured along the fluid congruence. The fluid congruence 
is usually inconvenient for computation, and so a second congruence utilizing a 
temporal gauge called the separable volume gauge [201 ^ s use d ( e -g-> the orthogonal 
frame for SH models). Let t and H be the proper time and the Hubble scalar 
associated with this congruence. The separable volume gauge allows a convenient 
computational time variable r, called the dynamical time, to be defined globally 



Perfect fluids with non-zero pressure have non-zero 4-acceleration when u a ^ n a . 
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Bianchi 
model 


Matter 


Extreme 
tilt sink 


Other sinks 
in set 


II 


4r < 7 < 2 


£ (II) 


No 


IV 


,6 < 7 < 2, J < £+ < 
5 ^ ! -< 7 <2, -i <S + <-i 


£_(IV) 
fL(IV) 


Yes for 7 < 4+ 3 J+ 

' 3 

No 


V 


^ <7<2 


M_ 


Yes for | < 7 < § 


VIo 


t<7<2 


f(VI ) 


No 


VI ft 


lh < 7 < 2 [*] 


[*] 


[*] 


VII 


|<7<2 


n(VIIo) 


No 


VUh 


5+lcT<7<2, -|<S + <0 
^-< 7 <2,-I<S + <-i 


£-(Vll,0 
f(VII ft ) 


Yes (V) for 7 < 
No 


VIII 


1 < 7 < 2 


Si (VIII) 


No 



Table 1. Summary of extremely tilted future attractors. [*] There 
are extremely tilted attractors in the Bianchi VI^ models, but not 
all of the sinks have been determined in this case to date; for 
— 1 < h < 0, h 7^ —1/9, we have that 7,1 < 6/5 (for h < —1, we 
have that % < 4/3). 



for the entire congruence of world-lines: 

< 2 > 

For an individual fluid world-line, along which the spatial coordinates are held 
fixed, we have the following relation between T and t: 



(3) dT = v 7 ! -v 2 dt . 

We can now express the proper time elapsed, AT, along this fluid world-line as an 
integral over r: 

(4) AT = T 2 - Ti = jT 2 v 7 ! - v 2 dt = J ' i ^1 - « 2 dr . 

Physically, r is directly related to the length scale £; H = and imply that 
i = £o(x l )e T . As I — > 00, so does r. In this paper we are interested in whether AT 
is finite or infinite as ti — + 00. If AT is finite and H —f 00 as t 2 — > 00, then the 
congruence is future incomplete: the fluid observers will reach infinite expansion 
within finite proper time. Therefore, in spite of the fact that these spacetimes are 
future geodesically complete |5T], the world-lines defined by the fluid congruences 
can approach null infinity so quickly that the proper time as measured by the fluid 
observer is finite. 

In this paper, we study perfect fluid cosmological models in which there are two 
naturally defined congruences, a geometrically defined geodesic congruence and a 
non-geodesic fluid congruence. In Appendix B we establish an appropriate set of 
boost formulae relating the physical variables in the two congruences. We study 
expanding spatially homogeneous tilted perfect fluid models, with an emphasis on 
future evolution with extreme tilt. In order to set the stage, we shall first discuss the 
Bianchi type V model in some detail. We shall then discuss spatially homogeneous 
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models with a cosmological constant and various models (Bianchi types VIIo, VH/t 
and VIII) in the absence of a cosmological constant. We then briefly discuss general 
spatially inhomogeneous models with a cosmological constant. We find that for 7 > 
4/3, expanding spatially homogeneous tilting universes come to a violent end after 
a finite proper time, reminiscent to that of cosmological models with a 'phantom 
energy' equation of state parameter which end in a big rip. We also find that for 
7 < 4/3) the tilt can become extreme at late times and give rise to an effective 
quintessential equation of state. 

To establish the potential connection with phantom cosmology and quintessence, 
we determine the effective equation of state in the models under consideration as 
seen by the fluid observer. In order to investigate the future asymptotic behaviour of 
the cosmological models we then determine the asymptotic behaviour of the tilting 
models in the fluid frame variables using the boost formulae derived in Appendix 
B. We also study the asymptotic dynamical behaviour of a class of LRS Bianchi 
type V solutions into the past, and show that they can be extended along the 
fluid congruence. We conclude with a discussion, with an emphasis on the physical 
interpretation of the models. 

2. Bianchi type V 

The irrotational Bianchi type V perfect fluid models were studied in JOj ; it was 
shown that the tilt can become extreme (in a finite time as measured along the 
fluid congruence). The global structure of Bianchi V models was studied in more 
detail in 0. It was shown that for models with 4/3 < 7 < 2, to the future, as T 
approaches some finite limiting value, the fluid world-lines become null (with respect 
to the geometric congruence), and the expansion, the shear, the acceleration, and 
the length scale I all diverge but the matter density (and curvature scalars) tend 
to zero; i.e., the models end at a 'conformal singularity' (this peculiar behaviour is 
depicted in Figs 7(ii) and 8 in 0; dynamically an observer is accelerating so quickly 
that their world-line becomes null near null infinity in finite proper time, signifying 
the onset of structural instability in the field equations in these models). 

This peculiarity occurs for ultra-radiative equations of state. Such an equation 
of state may not be physically appropriate at late times, although in j22| it was 
argued that the future is dynamically radiation dominated (but with an always 
sub-radiative equation of state). This perhaps indicates that models (or rather the 
governing equations) with an ultra-radiative equation of state may not be well- 
posed at late times (i.e., structurally unstable). However, we shall show later that 
a similar type of physical behaviour may occur for particular SH universes with a 
sub-radiative equation of state. 

2.1. Analysis. We shall first show that the fluid proper time is finite as the solution 
approaches the equilibrium point M~ for the case 4/3<7<2in the absence of 
a cosmological constant (see Fig. 10 in JO])- We shall make use of the evolution 
equations (denoted by HW) in JOj- There, the equilibrium point M~ is described 

by 

(5) £+ = E_=0, A = l, v = -l, 

with the density parameter O = from HW(2.13) and the deceleration parameter 
q = from HW(2.14). By linearizing equations HW(2.10) and HW(2.11) about 
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M , we heuristically obtain the decay rates as r — > oo for the Hubble scalar 
H = 9/3 and the quantity vl — v 2 : 

(6) H sa i/ e" T , v 7 !-^ 2 « (v 7 ! ~ " 2 )o exp 

which have been confirmed by our numerical simulations. Finally, it follows from 
Ijljl that the fluid proper time in the future asymptotic regime is 

(7) AT « ^ — exp ^-^^J d - - 

which is finite for 4/3 < 7 < 2, as claimed. This is consistent with the analysis of 
Collins and Ellis second table in Fig. 7] . 

2.2. Cosmological constant. Spatially homogeneous cosmological models with a 
positive cosmological constant were investigated using dynamical systems methods 
in extending the tilted LRS Bianchi type V analysis of JU] to the A ^ case. 
Unlike the case of orthogonal Bianchi models (including type IX), in which the de 
Sitter point is the only new equilibrium point, when extending these models from 
the A = case tilted models have additional equilibrium points, namely the de 
Sitter point with extreme tilt for 7 7^ 4/3, and with arbitrary tilt for 7 = 4/3. 

In the LRS case, the field equations and conservation equations give rise to a low 
dimensional constrained dynamical system when written in terms of dimensionless 
variables. Focusing on the future evolution of these models, all of the sinks were 
determined. The conventional de Sitter solution is a sink for 7 < 4/3. The only 
sinks in the case 7 > 4/3 are the de Sitter points with extreme tilt (v — ±1). 
The threshold value 7 = 4/3 is particularly interesting, as a non-zero cosmological 
constant then implies v' = for any value of v. In other words, for each value of 
the tilt variable v, there is a de Sitter equilibrium point |23| . 

Therefore, in general, for Bianchi type V models with 7 > 4/3, the tilt becomes 
extreme at late times. The fluid motion is no longer orthogonal to the surfaces 
of homogeneity, and the rest spaces of an observer comoving with the fluid need 
not be spatially homogeneous and the model will not seem to isotropize. (When 
following the normal congruence, on the other hand, the fluid will no longer look 
perfect.) 

It has been shown that expanding non-type-IX Bianchi models with a positive 
cosmological constant isotropize to the future (cosmic no hair theorem) |24| . This 
result applies to tilted models, but the isotropization is with respect to the con- 
gruence normal to the homogeneous symmetry surfaces - not the fluid congruence 
[2"5] . Thus in the Bianchi type V models the the space-time generically becomes de 
Sitter-like, in accordance with the cosmic no hair theorem, but since the tilt does 
not die away, isotropization of the cosmology does not occur with respect to the 
fluid congruence. It may be possible to extend beyond the apparent limiting point 
to a region with time-like symmetries, as discussed in [2]. 

3. Spatially Homogeneous Models 
In general, SH models are not asymptotically isotropic. 

3.1. Isotropic Asymptotia. The spatially homogeneous models isotropize to the 
future in the presence of a positive cosmological constant [21] • 
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3.1.1. Cosmological constant. Let us investigate a SH model with a cosmological 
constant and a tilted perfect fluid with p — (7 — l)/i. Consider the asymptotic 
behaviour of the following solution which corresponds to an equilibrium point of 
the system of governing equations: 

Solution: de Sitter (dS) with an extremely tilted perfect fluid. 
Asymptotic behaviour 26 : H = H , ^1 - v 2 « - v 2 ) exp (-^5^t) . 
This gives for AT: 

(*\ AT ~ rWIE^k ( 3 ^- 4 ^ A 

(8) J TQ —^ e ^{-— T ) dT ' 

which is finite for 4/3 < 7 < 2. This solution is the future attractor for 4/3 < 7 < 2 
and therefore represents the generic behaviour for this model. 

3.1.2. Intermediate behaviour. In the absence of a cosmological constant, in general 
SH models are not asymptotically isotropic. However, these models can spend a 
long time close to a Friedmann model corresponding to a saddle point :2 . Therefore, 
let us consider one non-tilting perfect fluid with p± — (j± — l)fi± and one tilting 
fluid with p = (7 — l)/x. We assume that the non-tilting fluid is dominant and has 
7_L < 7. Consider the equilibrium point: 

Solution: Friedmann F, with H w Ho exp(— §7±t) with extremely tilted second 
fluid. 

Asymptotic behaviour: \/l — v 2 ~ (VT— w 2 )o exp ^— ^--y 1 ^) • 
This gives for AT: 

In the class of SH models F is a saddle for 2/3 < 7^ < 2 and represents possi- 
ble intermediate behaviour. This means that for solutions spending a finite (but 
arbitrarily long) time close to the saddle F, exhibiting quasi-isotropic behaviour 
consistent with observations [5], the Hubble parameter as measured by the fluid, 
-Hfluid, can become arbitrarily large. Of particular interest are the cases when the 
non-tilting fluid is dust (71 = 1) or radiation (7^ = 4/3) for which AT is finite for 
7 > 14/9 and 7 > 8/5, respectively. 

3.2. Anisotropic Asymptotia. There will be different asymptotic end-states de- 
pending on Bianchi types (where the criteria may depend on 7). In all of the cases 
below we have one tilted perfect fluid. 

3.2.1. Bianchi type VIIo. The investigation of the type VIIo model is important for 
several reasons. First, the type VIIo model is the most general spatially homoge- 
neous model allowing for the flat Friedmann- Robertson- Walker (FRW) model as a 
special case. Both the Bianchi type I model and the VIIo generalise the flat FRW 
model; however, since the Bianchi type I perfect fluid model does not allow for tilt, 
only the Bianchi type VIIo model can be used to investigate the effect of tilt on 
the evolution of the universe close to flatness. In addition, the type VIIo model 
plays an important role in the Bianchi hierarchy. For example, the type VIIo is a 
special limit of both the 'most general' type VIII and type IX models. Also, unlike 
the other Bianchi models of solvable type, the type VIIo state space is unbounded. 
The type VIIq model is not asymptotically self-similar to the future because the 
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unbounded growth of one of the curvature variables leads to oscillations in the shear 
and the curvature variables [3J. 

The Bianchi type VIIo model with a tilted 7-law perfect fluid was studied in 
|18j using expansion-normalised variables and utilizing a formalism adapted to the 
time-like geodesies orthogonal to the hypersurfaces defined by the type VIIo group 
action. In particular, the future asymptotic states were determined. Regarding the 
tilt velocity, it was shown that for fluids with 7 < 4/3 (which includes the important 
case of dust, 7=1) the tilt velocity tends to zero at late times, while for a radiation 
fluid, 7 = 4/3, the fluid is tilted and its vorticity is dynamically significant at late 
times. For fluids stiffcr than radiation (7 > 4/3), the future asymptotic state is an 
extremely tilted spacetime with vorticity. 

For 4/3 < 7 < 2, the solution corresponding to the sink JSj has the asymptotic 
form 

(10) H^H e~ 2 \ Vl-v 2 ~ {Vl-v 2 ) Q exp 

This implies, 

which corresponds to the proper time of an observer from t = tq to r = 00 following 
the fluid congruence, there is a change of behaviour at 7 = 8/5. For models with 
7 < 8/5, we find that this integral diverges. However, for fluids with 7 > 8/5, we 
find that this integral is finite. Since asymptotically at late times a fluid with 7 > 
8/5 will be extremely tilted, this means that the fluid will reach future null infinity 
in finite proper time. This is consistent with the work of | 22| 2 . As discussed earlier, 
a similar phenomenon occurs for the LRS Bianchi type V model with 7 > 4/3 UJ. 

3.2.2. Bianchi type Vllh- For this model the appropriate solutions are vacuum 
plane waves with extreme tilt (see ^1 ^] for details). There are two possible 
asymptotic states depending on a parameter S + : 

(1) £-(VII h ): Stable for 6/(5 + 2S+) < 7 < 2, -1/4 < £+ < 0. 
Asymptotic behaviour: H w _ffoe~ < - 1 ~ 2S+ - )T , 

A 2 ~ , ft ^ ( 2(5 7 - 6 + 2 7 S + ) \ 

VI - V 1 PS (VI ~ t^JoCXp I — — T \ . 

This gives for AT: 

(2) S(VII h ): Stable for 3/(2 - £+) < 7 < 2, -1/2 < £+ < -1/4. 
Asymptotic behaviour: H ps H e~( 1 ~ 2S+ ) T , 

A 2~ f A aS ( 2(1 - 2S + )(2 7 - 3 - 7 S+) \ 

VI - , 2 « (VI - - 2 )o exp ^ 3|E + |(2- 7 ) T ) ' 

2 In particular, the dynamical future of radiation dominated Bianchi VIIo models, regarded 
as homogeneous perturbations of flat Friedmann model, was studied in 1221 . At late times it 
was found that v ~ £( 3 ~f~ 4 ), uj/H ~ ^( 9 T _10 ) ) signifying a growing vorticity instability. It was 
consequently found that the flat radiation filled Friedmann model is unstable to the development 
of vorticity. 
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This gives for AT: 



AT ps / ^^expf-(l^±j!^H dr. 



H "V 3|E+|(2- 7 ) 

In the case — 1/2 < E + < —1/4, the asymptotic state is actually a closed orbit, so 
in calculating the decay rates we have considered the appropriate averaged value 
as described in Section 4 of |16j . From the above expressions AT is finite whenever 
7 > (4-2E+)/3 for -1/4 < E+ < and 7 > 3/2 for -1/2 < E+ < -1/4. However, 
since for a given type VII/i model any —1/2 < E + < is possible, we have that: 
For tilting perfect fluid models of Bianchi type Vllh with equation of state parameter 
4/3 < 7 < 2, there exists a set of solutions of non-zero measure for which AT is 
finite. 

3.2.3. Bianchi type VIII. These models are asymptotically extremely tilted for 1 < 
7 < 2 and the asymptotic solution is an extremely Weyl-curvature dominated model 
(in the language of |27j). The solutions tend to a vacuum state but do so rather 
slowly: n/H 2 cx 1/r. We refer to f° r details. 



Asymptotic behaviour (1 < 7 < 2): H w H t* exp (— |r 



y/l - v 2 w {y/l- v 2 ) t 2 < 2 -^ exp (- 



2 1 

3(7-1) 
2-7 



This gives for AT: 

which is finite for 4/3 < 7 < 2. 

4. General Spatially Inhomogeneous Models with a cosmological 

constant 

It is known that for general spatially inhomogeneous perfect fluid models with 
a cosmological constant, the de Sitter solution with extreme tilt 3 is locally stable 
for 4/3 < 7 < 2 26] |^. The de Sitter asymptotic expansion is given in [2S1 - From 
Equations (3.43) and (3.28) of we have that 



v/l - ^« (vT^) exp f-h-± T ^ 



(13) H « Ho 

and thus 

which is finite for 4/3 < 7 < 2, as required. 

A comment on the physical ramification of this result is in order. We note that 
for 4/3 < 7 < 2, the generic behaviour is v — > 1, which implies that inflation does 
not isotropise an ultra-radiative fluid. Even if inflation is turned off after a certain 
number of e-foldings, -Hfl u id can become arbitrary large. However, as noted above, 
this result does not contradict the cosmic no- hair theorem |24j . since generically 
the space-time is dynamically de Sitter-like (in the normal frame). 



3 Here the tilt refers to the fluid tilt with respect to a congruence whose acceleration tends to 
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5. Effective equation of state 

We observe that the future dynamical evolution of the models under investigation 
(with 7 > 4/3) is reminiscent to that of cosmological models with a constant 
equation of state parameter y e g < 0, dubbed 'phantom energy', in which expanding 
universes end in a singularity called the big rip after a finite proper time |28| . To 
establish the connection with physics let us determine the effective equation of 
state, as measured by the fluid observer, in the models under consideration. 

For convenience, we first present the generic boost formula (see Appendix B) for 
the Hubble scalar: 

(15) ff fluid = TBH , 
where 

(16) B = 1 - \v a A a + \v a U a + Ir-^ofT) + |r- 1 fl a (T« a ) . 

A fluid version of J5J can be defined along an individual fluid world-line: 

l^7\ drfluid _ n 

K 1 ') ~dT~ ~ fluid ' 

Then the two dynamical times are related by 

(18) dT fluid = SdT. 

We note that a finite AT implies the divergence of -fffluid- Indeed, in the examples 
above -Hfluid diverges. 

Let us now consider that the fluid observer interprets the Universe as a fiat 
Friedmann universe with an effective equation of state. We define the effective 
density /i c ff and effective pressure p e s for the fluid frame by writing the Fried- 
mann equation as = 3i/| uid — fx e s, and the Raychaudhuri equation as eo-ff fluid = 
— -Hfluid - \ {Hcft + 3p c ff), as is standard. Equivalently, expressed in terms of Hubblc- 
normalized variables, we have fl e s = 1, <7fluid = ^(^+3P) e ff- We define the effective 
equation of state parameter 7 g as 

(19) 7cff = ^ + 1 , 



off 



so that 7eff is given by 



(20) 7off = f(g flui d + l). 

If the limit of Qfl u id as r — > oo exists, then it can be read of from the leading 
exponent of Hfl u id since 

(21) ff fl uid w (fffluid)oe- (9n » ,d+1)Tfluid , 

and Tfl U id = Bt in the limit. Thus if fffluid ~ (fffluid)oe fer , then gfl U i d + 1 = —k/B 
and 

(22) Te ff = "||. 

For the LRS Bianchi type V models 4 , on the approach to M~, we have that 

(23) k^J~ A \ B 1 



3(2-7)' 



4 The effective equation of state in the LRS Bianchi type V models can be calculated directly 



in fluid- frame variables (see Appendix 1X1. 
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and so 

(24) 7off = -(37-4) <0 for 7 >|; 

i.e., on approach to M~ the ultra-radiative fluid effectively behaves like a phantom 
energy in an isotropic and spatially flat spacetime 5 . 

Similarly, for the de Sitter asymptotic state we have that 

(25) * = ^i B 2 



2- 7 ' 3(2-7)' 
and coincidentally 

(26) 7cff = -(37-4) < for 7 >§. 

This is the general case in inhomogeneous cosmological models with a cosmological 
constant. As the de Sitter asymptotic state is approached the ultra-radiative perfect 
fluid effectively behaves like a phantom energy in an isotropic and spatially flat 
spacetime. Note that inflation does not stop the big rip from occurring. 

Let us now consider the SH models discussed above. For the Bianchi type VIII 
example, we have that 

(27) fc = 3 jf:- 4 ^ B = 1 



2(2-7) ' 2- 7 ' 
and again 

(28) 7o ff = -(37-4) < for 7 >§. 

We also note that for the Bianchi VIII model the future attractor is of extreme tilt 
for 7 > 1, and so < 7 e ff < 1 for 1 < 7 < |. 
For the Bianchi type VIIo example, we have 

(29) k=%=* t B = ' 



2-7' 3(2-7)' 
and 

(30) 7 cff = -(5 7 -8) < for 7 >§. 

We attribute this difference from the previous cases to the fact that Bianchi type 
VIIo models are geometrically more special than Bianchi type VIII models. For the 
first Bianchi type Vllh example, we have that 

(31) fc= 2(37-4 + 2S +)) B _4(l + S+ ) 



7 ' 3(2-7) 
and 



where E + satisfies — \ < S + < 0. And for the second Bianchi type Vllh example, 
we have that 

(l-2£+) 2 (2 7 -3) 



(33) k 



3|S+|(2- 7 ) 



^Note that for 7 = | , we obtain an effective equation of state corresponding to a cosmological 
constant. However, in this case the future attractors are not of extreme tilt and the asymptotic 
decay rates will differ. 
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and B is a complicated constant. But since B is always positive, we have that 
7cff < for 7 > §. Similarly, for the two-fluid example, we have 7 ff < for 

7 > 2 ^2+j^) ■ I n the generic Bianchi type VIII example the ultra-radiative perfect 
fluid effectively behaves like a phantom energy. In the more special examples, the 
requirement is that the equation of state must be equal to or higher than that of 
radiation. 

In order to compare the future asymptotic behaviour with the behaviour of 
various types of big rip singularities |281 2!) , in the next section we shall determine 
the asymptotic behaviour of the tilting models in the fluid frame variables. 



5.1. Quintessence. We also see from above that there are models which have 
extreme tilt asymptotically, but for 7 < | the effective equation of state is that 
of a cosmological constant (7 G ff = 0) or that of quintessence (0 < 7eff < §)• 
If 7 < |, then the future attractor may not be of extreme tilt, in which case the 
corresponding dynamics is very special. But, from Table 1, there do exist some 
extreme tilt sinks for 7 < |. For example, for Bianchi type VIII models there is a 
unique future attractor E\ which is of extreme tilt for all 7 > 1 , but from equation 
(|2"5)l we see that there are models with ^ < 7 < | that lead to < 7 c s < |, and 
thus generically mimic the dynamics of quintessence cosmological models. In the 
case of Bianchi models of type V, for | < 7 < | we again obtain quintessential 
behaviour (a similar behaviour will occur in type VIo and VI ^ models; however, in 
all of these cases the extreme tilt sink is not the only future attractor). For certain 
values of S + (see Table 1), the Bianchi type VII/j models also lead to quintessence. 
We note that spatially homogeneous and inhomogencous models with a cosmological 
constant cannot generically give rise to quintessential behaviour since there appear 
to be no future attractors with extreme tilt when 7 < | (although, of course, the 
models are already inflating). 



6. Fluid frame variables 

Let us next determine the asymptotic behaviour of the tilting models in the fluid 
frame variables using the boost formulae given in Appendix B. In particular, we 
shall present the de Sitter limit and the Bianchi type VIII future asymptotic limits 
as examples. 

For the de Sitter limit, the asymptotic form is given in (eqs. (3.22)-(3.28) 
and (3.43)). Note that the asymptotic form for v a in j2Ej, eqs. (3.27) and (3.28) 
can be simplified to 

(34) v a = vS + 0(e~ rT ) , l-v 2 =e- rt [(l-v 2 ) + O(e- rT )} , 

where r = 2(37 — 4)/(2 — 7). That is, the 0(e~ T ) term is zero. The outline of the 
proof is as follows; we choose a spatial frame such that v a — > (1, 0, 0) and we show 
that the 0(e~ T ) term is zero. The result is true for any Fermi-propagated frame. 
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With this simplification, we apply the boost formulae in the appendix to the 
asymptotic forms to obtain the fluid frame variables for the case 7 > | : 

(35) H &uid = ^—^^A^ + 0( e -4" + e^-^) 

(36) a£ uid = -ToJ^vSei" + [coe^' 1 * + 0(1 + e^' 1 ^) 

(37) ^ uid = ^Z^ro /I v ^ krT + [coef]^" 1 ^ + 0(1 + e^-^) 

all of which diverge as r — > 00, while 
(39) < d) < id , <d = 0(l + e M 



(38) < id = ^To^^e^ + [coer^-^ + 0(1 + e <^-Dn 



10 
7 • 



which diverge for 7 > 

The factors that contribute to the divergence (in the e^ rT term) are the TH 
term (affecting if fluid , ag uid and u% uid ), and the e terms (affecting H Ruid , tig uid 
and crg^. d ). We now look at the general situation for the existence of a singularity. 
We observe that, in order for a singularity to occur along the fluid congruence, the 
factor TH in Q must diverge sufficiently rapidly. The factor TH also appears in 
the boost formulae for H, a a and u a , so these variables for the fluid frame generally 
diverge, unless any cancellation occurs. This indicates that a singularity is generally 
characterized by the divergence of H, a a and u a . This leads us to speculate that 
a singularity along a congruence is primarily affected by the rate of expansion and 
the acceleration of the congruence rather than its shear and vorticity. Furthermore, 
the spatial curvature component a a plays a more significant role than n a /3 and the 
frame rotation fl a p plays a less significant role. 

We next present the asymptotic forms for the dominant variables for the Bianchi 
type VIII limit. From [T§], the asymptotic forms for the variables along the con- 
gruence normal to the SH hypersurfaces are 

(40) HttH Q e~* T , a+ « ±H Q e~i T , n+ » (n+) r"i , 

(41) 

1-t) 2 w (1 -w 2 ) r^e" 1 fer iT , V! fa 1 + 0(1 -v 2 ) , v 2 ,v 3 w 0(r3e-i T ) 

(42) t7 12 ,a 13 «0(r-3e- 3T ) , CTl = (cr 2 + a 2 + n 2 _ + n 2 )tf 2 = ©(r" 1 ) , 

as t — ► 00 , valid for 1 < 7 < 2. Applying the boost formulae in the appendix, we 
obtain 

(43) iJfluid « -^—T Q H Q T-^)e^r T 



(2 - 7)i/fl uid , Wfl uid ~ (27 - l)-fffl u id , (cr + )fi uid w (37 - 4)if flt 



(44) 

a fluid 

(45) (n+)fluid » ro^+jor^e^r 1 , 

as t — ► 00, while other fluid frame variables are less dominant. 

As in the de Sitter example, for | < 7 < 2, the factor TH contributes to the 
divergence of £/fl u id, flfl U id and Wg uid , while the e terms contribute to the divergence 
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of H fl u id, u nuid and (cr+ )fl u id- In addition, the diverging factor Tn+ contributes to 
the diverging (ro+) n uid, and Ta+ contributes to (o+)fl u id- 

From the asymptotic behaviour of the extremely tilting models in the fluid frame 
variables we have found that the length scale and the Hubble parameter (i.e., the 
expansion) diverge in finite proper time, and the singularity is consequently a future 
sudden singularity in the sense of Barrow |29| . In addition, the shear and the 
acceleration of the fluid congruence also diverge (while the energy density tends 
to zero). That is, the fluid congruence experiences a 'finite kinematic' singularity, 
and shares many of the physical properties of a big rip singularity ,28 . A 'finite 
kinematic' singularity is the analogue of a 'conformaP singularity [5J |§] or a 'finite 
density' singularity [2j (in which the energy density tends to a finite value) that 
occurs in the Bianchi type V models. In particular, at a 'finite kinematic' singularity 
all of the Weyl invariants remain bounded (although components of the Weyl tensor 
can diverge). 



As noted earlier, a similar situation can occur to the past. As a final example, 
let us consider the asymptotic dynamical behaviour of a class of Bianchi type V 
solutions. Collins and Ellis |3 showed that there exists a set of typical LRS Bianchi 
V solutions that, into the past, can be extended along the fluid congruence, from a 
spatially homogeneous region of spacctimc into a stationary region (see their figures 
3-7). Along the normal congruence, however, an initial singularity is reached, with 
the limits 



and H m ir e~ (1 ~ 2S + o)r , T « r e~ (1-2S + a)T . The proper time elapsed along the 
fluid congruence up to the Cauchy horizon is 



7. LRS Bianchi type V extendible solution 



(46) 



(E+o.l + £+o,l), -K£+o<0, 



(47) 



AT = 



TO X 




1 



e 2(l-2S +0 )r dr 



— OO 



TH 



— oo 



which is finite. 

Here, we apply the boost formulae to the asymptotic forms to check that the fluid 
frame variables are bounded at the Cauchy horizon, consistent with the extendible 
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solution: 6 
(48) 



(2 - 7 )S+2 + 2S +0 - 2 7 H , ^ /4 i_ aE+o)7 



Hfiuid = -± " ™ - + e 

37(2 - 7) r 

(49) 

(aiWd - |(1 + S+o)(l + S+o - ^+o)^ + 0(e^- 2S +^) 

' 1 

(50) 

(7-l)[(2-7)S + g + 2E +0 -2 7 ] g o(e(1 _ 2S+0)n 

^ijfluad- 7(2-7) r + 1 ' 

(51) 

, x (3 7 - 4)(2 - 7 )S + g - (6 7 2 - I87 + 8)S +0 + 7(37 - 2) H (1 - 3S+0 ) n 

(^+)fluid = ^— ^ : 7r + U ( e ) i 

67(2-7) r 

as t — > — 00, while other components of the fluid frame variables are identically 
zero. 

Therefore, in this case, although the proper time as measured by an observer 
is finite, there is a congruence along which the evolution can be continued into 
the past. This raises the question of whether, for any particular exact cosmological 
model or corresponding asymptotic state, there exists a congruence whose observers 
cross a Cauchy horizon into a space-time region not accessible by the observers along 
the standard congruence of the model considered. In particular, it is of interest 
to determine whether there is a congruence which can be continued beyond the 
asymptotic state in any of the examples discussed above in which an extreme tilt 
future asymptotic state is approached in finite proper time. 

8. Discussion 

We have found that for 7 > | a 'finite kinematic' singularity develops in which 
the nature of geometrically defined surfaces change from spacelike to null in a finite 
time as measured by the fluid. The field equations are structurally unstable, and the 
question consequently arises as of whether these models are mathematically well- 
defined. In any case, to fully understand the behaviour of these models and their 
physical properties, the dynamics need to be studied using a formulation adapted 
to the fluid (i.e., utilizing a fluid-comoving frame). 

On the other hand, this mathematical behaviour might lead to some interesting 
physics. Expanding universes that come to a violent end after a finite proper 
time have arisen in a different context. Models with a constant equation of state 
parameter 7 < 0, dubbed 'phantom energy', lead to a singularity called the big rip 
|28| . In this paradigm, during the cosmic evolution the length scale grows more 
rapidly than the Hubble distance and consequently blows up in a finite proper 
time, and is typically characterized by a divergent pressure and acceleration. As 



6 The asymptotic forms for and A need to be obtained explicitly up to the order 

e 2(i-2£ + o)-r^ anf j p U p £ Q £ ne or( j er e (i-2S +0 )T^ as cance ]l a tions occur at leading orders to give 
the eventual bounded limits. Heuristically, we obtain the asymptotic forms as follows: the ansatz 
£+ = S+o + S+ie fc i T + A = 1 + S+o + Aie k z T + is substituted into the £+- and A- 
evolution equations and the constraint equations to give k\ = k^ = 2(1 — 2E_|_o), and and A\ 
explicitly in terms of S+o and To. This ansatz is then substituted into the T-evolution equation 
to give T = r e~ ( - 1 ~ 2s +°) T + rie( 1_2E +°) T , where Ti is given explicitly in terms of S +0 and T . 
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the big rip singularity is approached, both the strong and weak energy conditions 
are violated. 7 The details of the pathological behaviours depend on the particular 
phantom cosmological model under consideration, and less violent singularities such 
as future sudden singularities are possible 29 . 

To establish the connection with phantom cosmology we determined the effec- 
tive equation of state in the models under consideration. We found that as the 
asymptotic state is approached the ultra-radiative perfect fluid effectively behaves 
like a phantom energy in an isotropic and spatially flat spacetime. It is important 
to note that the energy conditions of the perfect fluid are nowhere violated. In 
order to compare the future asymptotic behaviour with the behaviour of various 
types of big rip singularities ,25j and future sudden singularities (2H| we computed 
the future asymptotic dynamical behaviour of the tilting models. 

Let us discuss the physical consequences of this dynamical behaviour in a little 
more detail. The models presumably spend a period of time approximately isotropic 
(i.e., close to a flat Friedmann saddle point), with a very small (but non- vanishing) 
tilt. Thereafter, the models begin to evolve away from isotropy. Since the tilt is 
non-zero, for 7 > § the models generically evolve towards an asymptotic state with 
extreme tilt. As the tilt becomes extreme, observers moving with the tilting fluid 
will experience dynamical behaviour mimicking that of a phantom cosmology. We 
note that such behaviour seems to be consistent with galaxy, CMB and supernovae 
observational data |30| . Moreover, unlike in conventional phantom cosmology, in 
the models studied here there is no need for any exotic forms of matter; conventional 
matter which is tilting suffices. 8 Indeed, other pathologies, such as the existence 
of ghosts, are avoided in the models described here. 9 In addition, as noted above, 
due to the existence of future attractors with extreme tilt the dynamical behaviour 
described here is generic. 

We have also found that there are tilting SH perfect fluid models with 7 < § 
which have extreme tilt asymptotically and have an effective quintessence equation 
of state (0 < 7efj < |) EH- Unlike the case of the cosmologies with 7 > § discussed 
above in which the phantom behaviour is generic, the quintessential behaviour is 
typical (i.e., there are other future asymptotic behaviours possible) and model de- 
pendent (i.e., it does not occur in all Bianchi models). The details of the extreme 
tilt sinks in different Bianchi models for 7 < | is given in Table I, and the corre- 
sponding ranges of 7 that lead to < j c s < | and thus mimic the dynamics of 
quintessence cosmological models is discussed in Section 5. 

In future work we shall further study the physics of the models with tilt, from 
the perspective of the observers moving with the fluid matter congruence, using 
the formalism developed in this paper. We note that the tilting models studied 
here agree with observations in the same sense that the phantom cosmologies and 
quintessence models agree with observations, and thus deserve further scrutiny. In 
particular, we will investigate the effect of tilt on cosmic microwave background 



7 We are only considering perfect fluid models here; it may be of interest to investigate more 
realistic fluids with viscosity. 

8 In a braneworld approach, accelerating universes can also result without a cosmological con- 
stant or other form of dark energy 1321 . 

9 This is also the case in alternative models to phantom cosmology which result from alternative 
theories of gravity, theories with non-minimal couplings, and models in which the dark energy and 
quintessence field interact )33| . 
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radiation observations. It is possible that these tilting models may offer an ex- 
planation for the low quadrupole amplitude and other anomalies on large angular 
scales 31 found in the WMAP data [SB] . 
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Appendix A. LRS Bianchi type V models in fluid-frame variables 

In the LRS Bianchi type V case, the evolution equations in the fluid frame can 
be presented as a system of two ODEs, simplified from the equations for general 
spatially inhomogeneous models 10 by setting 

(52) di=vd Q , d 2 =d 3 = 7 £ Q/3 =diag(-2£ + , £+,£+), 

(53) A Q/3 = 0, A a =(A,0,0), r a = (r,0,0), U a = (U, 0, 0) , 

(54) R a =W a =0, P = 0, Q a — , n Q/3 = . 
The resulting system is 

(55) v' = (2E+ + 3 7 - A)v 

£'+ = -3(1 + S + )G: 1 [G+S+ + 2( 7 - 1)V] 

—^[(9 + v 2 )Y? + + 2(6 7 - 5)i> 2 £+ + (3 7 - 2)v 2 } , 
3vG_ 



(56) 




where 




(57) 


A 


(58) 


a 


(59) 


b 


(60) 


c 


(61) 


G± 



-b - \/b 2 - Aac 
2a 

-|(37-2) 

[-3G_ + (2 + G_)(1 + £+)]/?; 

2 



(2- 7 )(l + E + ) 2 + 3( 7 -l)G_(l + E + ) 



The other fluid-frame variables are given by 

(62) {7 = 3(7-1)?; 

(63) - = q+ 1 = -^-[G + E + + 2(7 - 1)V] - -^H(9 + v 2 )Z + + (3 7 - 2)v 2 ] 
v G- 3v&- 

(64) K = 2(7 - 1)(1 + Z+)v 2 + A 2 + |u(l + T, + )A 

(65) Q = 1-'E 2 + -K. 

On approach to the M~ point, v tends to —1 and the fluid-frame variables are: 

(66) £+ = -±(3 7 -4), A= f(2- 7 ), 17 = -3( 7 - 1) , r = |(3 7 -4), 
which lead to 

(67) fi = P = 0, £ 2 = ±(3 7 -4) 2 , K = |(3 7 -2)(2- 7 ) , o = -I(9 7 -10). 



10 See I ' j Ci I for the general equations in component form. 
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Many of the calculations are more straightforward in fluid frame variables. For 
example, using these equations we can easily calculate the asymptotic decay rates 
(10) and (11) on the approach to M~ . We can also calculate the effective equation 
of state (assuming that the effective solution is spatially flat and isotropic, so that 
K cS = = S off , O e ff = 1 and q Ruid = |(fi + 3P) off ): 

(68) 7off = -(3 7 -4) 
which is negative for 7 > | on approach to M~. 

Appendix B. Boost Formulae 
The commutation functions are given by: 

(69) [e a ,e 6 ] = 7 c afc e c 
which can be decomposed in the following way: 

(70) 7 ° 0a = ii a 

(71) 7 V = -H8j - <rj - - SI J 

(72) 7% = 2u; a/3 

(73) ^ af3 = e^vT + a a S^ - apS^ 

where a a p is symmetric and traceless; n a p is symmetric; u> a p and Q a p are 
antisymmetric. For convenience, we also give the expressions for the compo- 
nents in terms of the commutation functions (angled brackets denote trace- 
free symmetrization): 

(74) H = -jT'V ' a a/3 = 7(/30a> , + ^af3 = 7[/30a] 

(75) n^ = ^Y) p] a a = -±-f lia . 

The boost formulae for the orthonormal frame vector fields are: 

(76) e = re + 

(77) e a = Tv a e + £ Q % 
where 



(78) B a » 



s a " + Yqr[ VaVl1 



v 2 = 



The corresponding boost formulae for the commutation functions are: 
7° 0Q = T[e (Tv a ) - e Q (T)] - Ty^B^) - e a (Tv»)} 

+ (TB a » - r 2 v a v»)(T 7 ° 0f , - r^ 7 <v) 

(79) + Tv^B a v (Tj% u - Tvrf^) 

fV = -r^[eo(rt; a ) - e a (r)] + B/[e (B a ») - e a (Fv^)} 

+ (TB a » - rV^X-r^Vcv + iWV) 

(80) + WKAV + 
7% = 2Te [a (Tv m ) - 2rv^e [a (Bpf) 
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+ 2r V[a £ / /(r 7 ( V - rivyV) 

(81) +B^B p v {Y 1 % v -Tvtfuv) 
j\p = -2r^e [Q (r^]) + 2B^e [a (Bpf) 

+ 2r V[Q s / /(-rv'7 ( V + ^Vo,*) 

(82) + B^Bp^-Tv^^ + BJ 7 V) 

The corresponding boost formulae for the components of the commutation 
functions are then: 



(83) H 



r 4 

- rV(w„ <a + n^)^ + 2rV^ (aV/3) 

r 3 
r 3 

(84) + r 2 ?/^^) - Yip^^v^p) 

r 2r 

+ = e^Pi^]) - ^-^v [/3 e a] (r) + ^-jV { pe (Tv a] ) 

T 4 v 2 

+ r(u; a/3 + n a(3 ) + — _ _ + n^ju^ + 2Y 2 v tx u lx[a v l3] 
r 3 

(85) + r 2 a [Q v /3] + V 2 u [a v l3] 

r 

w a /3 = e^^j) - ^pY«[/3e a ](r) 

- r 2 ^ M 0- M [a«/3] 

2r 3 

r 3 

(86) - T 2 a [a v p] - r 2 u [aVf3] 

r r 

"a/3 = p^T e ^(a|^e /1 (rv| / j ) ) - -—e^^je^rv,,) 
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r 3 
2r 3 

r 4 

(87) - rVe^Wfla" - T 2 e Ma v ^ 

r r 

r 3 

- 2THv a + rV<7 a „ - ^—^j IJa ,vi i v v v a 

+ rVo^ + n Q „) - 2rV^ M 

p4 w 2 

+ (2 + rV)a a - — - a^t; a 

(88) + rV«a - I^i^wa 

r 

u a = e (rw a ) - ^-pY?; Q eo(r) 



+ + rVff^ - rr——<T^ v V^V U V ct 



r + r 

+ rV(uv + n aM ) - 2rV^ aM 

+ r Ve^n**^ + rVa Q - I%tA; a 

r 3 

(89) + r 2 u a - _^^ Q 

The energy-momentum tensor T a ft has a standard decomposition relative 
to a congruence of world-lines: 

(90) Too = fj, , T 0a = -g a , T a/3 = p5 a/3 + 7r a/3 , 

where \x is the total energy density, q a is the energy flux, p is the isotropic 
pressure and -K a p is the trace-free anisotropic pressure, all relative to the 
above congruence. These variables are called the source terms. 

Boosting the orthonormal frame results in the following boost formulae 
for the source terms: 

(91) A = H + r V (ji+p)- 21%^ + T 2 tt^v u 

(92) p = P + ±r V (n+ P )- fr 2 ^ + irV^t/y 

r 2 (2r + 1) 

g« = -T 2 {ji + p)v a + rg a h — — q v ,v ' J ' Va 
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(93) 



r 3 



a 



n a (3 = T 2 (fi + p)v {a v p) - 2Tv {a q p) 



2r 3 



r + i 

r 4 



(94) 



2r 2 



(r + i) 2 



2 f- v 



References 



1. G.F.R. Ellis and M.A.H. MacCallum, Comm. Math. Phys. 12 (1969) 108. 

2. Dynamical Systems in Cosmology, eds: J. Wainwright and G.F.R. Ellis, Cambridge University 
Press (1997). 

3. A. A. Coley, Dynamical Systems and Cosmology, Kluwer, Academic Publishers (2003). 

4. J.D. Barrow and D.H. Sonoda, Phys. Reports 139 (1986) 1. 

5. A.R. King and G.F.R. Ellis, Comm. Math. Phys. 31 (1973) 209; G.F.R. Ellis and A.R. King, 
Comm. Math. Phys. 38 (1974) 119. 

6. C.G. Hewitt, R. Bridson, J. Wainwright, Gen. Rel. Grav. 33 (2001) 65. 

7. I.S. Shikin, Sov. Phys. JETP 41 (1976) 794. 

8. C.B. Collins, Comm. Math. Phys. 39 (1974) 131. 

9. C.B. Collins and G.F.R. Ellis, Phys. Rep. 56 (1979) 65. 

10. C.G. Hewitt and J. Wainwright, Phys. Rev. D46 (1992) 4242. 

11. D. Harnett, Tilted Bianchi type V cosmologies with vorticity, Master's thesis, University of 
Waterloo (1996). 

12. S. Hervik, Class. Quantum Grav. 21 (2004) 2301. 

13. A. A. Coley and S. Hervik, Class. Quant. Grav. 21 (2004) 4193. 

14. J.D. Barrow and S. Hervik, Class. Quant. Grav. 20 (2003) 2841. 

15. A. A. Coley and S. Hervik, Class. Quant. Grav. 22 (2005) 579. 

16. S. Hervik, R.J. van den Hoogen and A. A. C oley, Class. Qua nt. Grav. 22 (2005) 607. 

17. W.C. Lim, R. J. Deeley and J. Wainwright, gr-qc/0601040 . 

18. S. Hervik, R.J. van den Hoogen, W.C. Lim and A. A. Coley, Class. Quant. Grav. 23 (2006) 
845 gr-qc/0509031 . ^ 

19. S. Hervik and W.C. Lim, to appear in Class. Quant. Grav. (2006) gr-qc/0512070 . 

20. C. Uggla, H. van Elst, J. Wainwright and G. F. R. Ellis, Phys. Rev. D 68 (2004) 103502. 

21. A.D. Rendall, Math. Proc. Camb. Phil. Soc. 118 (1995) 511. 

22. J.D. Barrow and F.J. Tipler, Nature 276 (1978) 453. 

23. M. Goliath and G. F. R. Ellis, Phys. Rev. D 60 (1999) 023502. 

24. R.M. Wald, Phys. Rev. D28 (1983) 2118. 

25. A.K. Raychaudhuri and B. Modak, Class. Quant. Grav. 5 (1988) 225. 

26. W.C. Lim, H. van Elst, C. Uggla and J. Wainwright, Phys. Rev. D 69 (2004) 103507. 

27. J.D. Barrow and S. Hervik, Class. Quant. Grav. 19 (2002) 5173. 

28. R.R. Caldwell, Phys. Lett. B 545 (2002) 23; B. Mclnnes, JHEP 0208 (2002) 029; R.R. Cald- 
well, M. Kamionkowski, N.N. Weinberg, Phys. Rev. Lett. 91 (2003) 071301; L.P. Chimento 
and R. Lazkoz, Phys. Rev. Lett. 91 (2003) 211301: A. Coley, S. Hervik and J. Latta, 
|astro-ph/0503169"| . 

29. J.D. Barrow, Class. Quant. Grav. 21 (2004) L79; J.D. Barrow and C. Tsagas, Class. Quant. 
Grav. 22 (2005) 1563. 

30. A.G. Ricss et al., Ap. J. 607 (2004) 665; U. Seljak et al., Phys. Rev. D 71 (2005) 103515. 

31. P.J.E. Peebles and B. Ratra, Ap. J. 325 (1988) L17; C. Wetterich, Nuc. Phys. B 302 (1988) 
668 ; R.R. Caldwell, R. Dave, and P.J. Steinhardt, Phys. Rev. Lett. 80 (1998) 1582. 

32. G. Dvali, G. Gabadadze and M. Porrati Phys. Lett. B 485 (2000) 208; V. Sahni and V. Yu 
Shtanov, J CAP 11 (2003) 014. 

33. L. Amendola, Phys. Rev. D 62 (2000) 04 3511; G. Huey and B.D. Wandelt, astro-ph/0407196 ; 
S. Das, P.S. Corasanuti and J. Khoury, astro-ph/0510628 . 

34. T. R. Jaffe, A. J. Banday, H. K. Eriksen, K. M. Gorski and F. K. Hansen, Astrophys. J. 
629 (2005) LI |astrc-ph/0503213|; T. R. Jaffe , S. Hervik, A. J. Banday and K. M. Gorski, 
submitted to Astrophys. J. astro-ph/0512433 . 



FLUID OBSERVERS AND TILTING COSMOLOGY 



21 



35. A. de Oliveira-Costa, M. Tegmark, M. Zaldarriaga and A. Hamilton, Phys. Rev. D 69 (2004) 
063516 astro-ph/0307282 . 

36. H. van Elst, manuscript www.maths.qmul.ac.uk/ hve/onf compo-h.pdf 

Department of Mathematics & Statistics, Dalhousie University, Halifax, Nova Sco- 
tia, Canada B3H 3J5 

E-mail address: herviksOmathstat .dal . ca, wclimOmathstat . dal . ca, 
aac9mathstat .dal . ca 



